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ABSTRACT 

The surface density of submillimeter (sub-mm) galaxies as a function of flux, usually termed the source 
number counts, constrains models of the evolution of the density and luminosity of starburst galaxies. At the 
faint end of the distribution, direct detection and counting of galaxies are not possible. However, gravitational 
lensing by clusters of galaxies allows detection of sources which would otherwise be too dim to study. We 
have used the largest catalog of sub-mm-selected sources along the line of sight to galaxy clusters to estimate 
the faint end of the 850 /im number counts; integrating lo S = O.lOmJy the equivalent flux density at 850/im 
is vli, = 0.24 ± 0.03 nW m"^ sr"' . This provides a lower limit to the extragalactic far-infrared background and 
is consistent with direct estimates of the full intensity from the FIRAS. The results presented here can help to 
guide strategies for upcoming surveys carried out with single dish sub-mm instruments. 
Subject headings: cosmic background radiation - gravitational lensing: strong - submillimeter: galaxies 
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1. INTRODUCTION 

Measurements of the sources which comprise the cos- 
mic infrared background (CIB) radiation provide constraints 
on the most vigorous epoch of star formation in the uni- 
verse. Deep submillimeter (sub-mm) surveys using the Sub- 
mm C ommon-User Bolometer Array (SCUBA: [Holland et alJ 
11999^ have detected a population of high redshift, ultra- 
luminous sub-mm galaxies with star formation rates ap 



2 mm (for example 


Lagache et alJl2003t iNegreUo et alJl2007l: 


IValiante et alj|2009 


). The counts observed at 850/im play a 



proac h ing lOOOMfD per yea r (fSmail et al.|ll997t iHughes et al 

lBlainetal.1 11991 IBm^ 



19981 iBarger et all fT998t 



ger et al 



19991 lEales et a l. 2000; Cowie et a l. 2002; Scott et al.l | 2002 : 



Chapman etal. 2002; Smail et al. 2002; Webb et al. 2003 
Borysetal. 2003; Coppinetal. 2006; Knudsen et al. 2006 
Knudsen et al.l 12008). Such surveys have been very success- 
ful in characterizing the number counts of 850/im sources 
with fluxes > 2mJy, and have resolved as much as 50 % 
of the extragalactic background light at these wavelengths 
(e.g. [Coppin et al. 2006). However, due to SCUBA's confu- 
sion noise of ^ 0.5 mJy at lu and a source density which 
rises steeply with decreasing flux, obtaining constraints on 
the number counts of sources at the low flux es correspon ding 
to the bulk of the CIB has been difficult. BlainI (1199 7^ ad- 
vocated using gravitational lensing by galaxy clusters to am- 
plify the brightness of these sources above the James Clerk 
Maxwell Telescope (JCMT) confusion limit. This method 
was first u s ed in pioneering work with SCU BA data by 
Smail et al.' (1997*) and subsequently by others (Sm ail et"an 
2002; Cowie et al. 2002; Knudsen et al. 2006; Knudse n et aP 
2008), to measure the background contribution from faint 
point sources at high redshifts, as well as the shape of the 
point source counts down to ^ 0.1 mJy. Unfortunately, these 
studies suffer from large Poisson noise due to a small num- 
ber of sources in their sample as well as challenging analysis 
issues; a larger data set is desirable. 

Models of the evolution of star-forming galaxies are fit 
to counts of galaxies at wavelengths ranging from 24 /im to 
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particularly important role because they span a broad range in 
redshift, across which episodic starburst activity is believed to 
have varied widely. It is therefore particularly important that 
the counts be extended to faint enough fluxes to capture the 
bulk of the CIB, as has now been done at shorter wavele ngths 
by BLAST (iPatanchon et alJl2009HMarsden et al.ll2009l) . 

The slope of the 850/im source surface density must 
break to a fl atter slope than is measured for bright sources 
(e.g. Copp in et alJ |2006) in order that the total flux re- 
main consistent with absolute mea surements of the CIB 
jPuget et al.lll'9961: iFixsen et alj|1998h . The precise shape of 
this break will constrain models of the evolution of num- 
ber density and luminosity of fainter infrared-luminous dusty 
galaxies and will also provide clues about the relation of this 
population to more ordinary, luminous sub-mm galaxies. 

In this paper, we utilize the SCUBA Galaxy Cluster Survey 
from'Zemc ov et all (l2007h . the largest 850/im galaxy cluster 
data set to date, and a model of the lensing in each cluster to 
form a catalog of sources on the line of sight through the clus- 
ters. These sources allow the most reliable constraints avail- 
able on the faint source contribution to the CIB at 850/im. 
Knowledge of the amplitude and shape of the source density 
spectrum will help set the stage for new sub-mm experiments 
like SCUBA-2, LABOCA, and HERSCHEL which are capable 
of mapping hundreds of lensing clusters over unprecedented 
sky areas. 

2. DATA SAMPLE AND ANALYSIS PIPELINE 
2.1. Su rvey sources 

More than 40 clusters were mapped over SCUBA's oper- 
ational lifetime with integration times ranging from Iks to 
> lOOks; all those archival SCUBA data taken as part of 
programs to rnap galaxy clusters comprise the data set pre- 
sented in IZemcov et alj (12007). The initial catalog used in 
thi s analysis con s ists of the 850 /im sources listed in Table 5 
of IZemcov et alJ (l2007h : that paper gives a complete discus- 
sion of the data sample, low level analysis, map making and 
source extraction techniques. 

It is thought that most blank field ex tragalactic SCUBA 
sources are at high redshift (iLagache et al.i,2005) and that their 
flux is thermal dust emission heated by an evolving combi- 
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TABLE 1 

Cluster Fields, Redshifts and Velocity Pispersion 
Parameters used in this Analysis. 



Cluster 


z 


Vz 


Reference 


Abell 209 


0.209 




Mercurio et al. (2003') 


Abell 370 


0.375 


1364;50 


Struble et al. (1991) 


Abell 383 


0.187 


900+ji; 


Smith et al. (2005) 


Abell 478 


0.088 




Wu et al. (1999) 


Abell 496 


0.033 




Faddaetal. (1996) 


Abell 520 


0.199 


1250:|«^ 


Proust et al. (2000) 


Abell 586 


0.171 


1 161+'"' 


CvDriano et al. (2005) 


Abell 773 


0.217 




Smith et al. (2005) 


Abell 851 


0.407 


i08i;3[ 


Goto et al. (2003) 


Abell 963 


0.206 




Smith et al. (2005_) 


Abell 1689 


0.183 


1990+100 


Lokas et al. (200^ 


Abell 1835 


0.253 


1210+*" 


Smith et al. (2005_) 


Abell 2204 


0.152 


1029!^2 


Pimbblet et al. (2006) 


Abell 2218 


0.176 


1070+1 


Smith et al. (2005) 


Abell 2219 


0.226 


902+\° 


Smith et al. (2005) 


Abell 2390 


0.228 


iiooi« 


Natarajan et al. (2004) 


Abell 2597 


0.085 


776!1»1 


CvDriano et al. (2004) 


C10016+16 


0.541 


984+130 


Poggianti et al. (2006) 


CI J0023+0423A 


0.827 


41 S+102 


"Lubinetal. (1998) 


C10024.0+1652 


0.390 


1000+™ 


Natarajan et al. (2004) 


CIG J0848+4453 


1.27 




Rosati et al. (1999) 


MS 1455.0+2232 


0.258 


1032+if 


Borgani et al. (1999) 


CI J 1604+4304 


0.895 


982!io« 


Crawford et al. (2006) 


CI 2244-0221 


0.330 


600i« 


Natarajan et al. (2004) 


MS 0440.5+0204 


0.190 




Borgani et al. (1999) 


MS 045 1.6-0305 


0.550 


1330!iii 


Borgani et al. (1999) 


MS 1054.4-0321 


0.823 




Tranetal. (1999) 


RX J1347.5-1145 


0.451 


1500!160 


Fischer etal. (1997) 



nation of star formation and active galactic nuclei. We ex- 
pect this holds for the sources in this work with two classes 
of exception: those sources known to be situated inside the 
cluster under investigation (e.g. Abell 780-1) and those ap- 
parent sources which may instead be images of the Sunyaev- 
Zel'dovich (SZ) effect (e.g. Abell 1689-2), which is positive 
at 850 /im. The former are radio bright so cross identifica- 
tions are easily found in the literature; such sources are read- 
ily separated from the high-redshift sample. The latter class 
of source is necessarily near the center of the cluster where 
background sources are demagnified. We expect no sources 
behind the cluster center to be detected and treat all central 
sources as cluster members. These two classes of sources are 
discussed in detail in Zemcov et al. (2007). In this work, we 
assume that both types are sources in the cluster (that is, not 
gravitationally lensed by it). The flux from such sources can 
simply be summed without magnification to yield their con- 
tribution to the total flux density . 

The full Zemcov et alJ (|2007|) sample contains some maps 
which are extremely noisy. We excise from the full sample 
all maps where the pixel rms in a SCUBA beam is > lOmJy 
and those maps where the data are obviously corrupted by 
instrumental problems. We are left with a subsample of 28 
clusters with which to measure the contribution to the CIB; 
these are listed in Table[T] 

A direct estimate of the total background arising from 
lensed sources can be made without the intermediate step 
of constructing demagnified number counts by exploiting the 
fact that gravitational lenses are like any other optical system 
in that the etendue of the system, AVl, is conserved; this means 
that, integrating over all f2, the total flux is equal before and 
after passing through the lens. The flux density can therefore 



be computed by integrating all of the flux in the sources in 
this survey and dividing by the area observed; this procedure 
yields the CIB lower limit vl^, > 0.19nW m"^ sr"'. Unfor- 
tunately, this method is limited because the and depth ob- 
served by SCUBA in these 28 clusters is significantly smaller 
and shallower than that required to contain all the flux from 
the reported FIRAS background. We therefore expect this con- 
straint to be smaller than the result from a more comprehen- 
sive analysis. Furthermore, it is difficult to assign an uncer- 
tainty to this approach as the survey completeness is not a sim- 
ple function of the measurement error in the image plane and 
certainly depends on the true shape of the source counts. We 
therefore employ a gravitational lens model to make progress. 

2.2. Gravitational Lensing Demagnification 

We model the lensing of each cluster using the mass dis- 
tributio n associated with t he Navarro-Frenk- White (NEW) 
model dNavarro et alJI 19961) since that has been shown to be 
a very good description of systems whose mass is dominated 
by dissipationl ess dark matter. Th e NEW model is discussed 
in wo rks like IWrig ht & Brainer dI (2000) and Li & Ostrike^ 
( I2OO2I) which provide all of the materials necessary to cal- 
culate the deflections and magnifications of idealized lensing 
clusters as a function of the concentration parameter c{z) and 
velocity dispersion (Ty. 

The background source counts are determined from the cat- 
alog as follows. Eirst, the catalog fluxes and the survey com- 
pleteness are corrected for the effects of gravitational lensing 
using the lens model. The de-lensed flux counts are then com- 
pleteness corrected and summed in bins to determine the mea- 
sured cumulative source counts A^meas(> S). Simulations are 
performed using different input source count models to com- 
pute the effects of flux boosting (Eddington bias; see, e.g., 
lEales et ani2000h on the true, underlying background source 
counts A'true(> S); comparison of the simulated, flux-boosted 
counts A^sim(> S) with Nms:as(> S) yields the corresponding 
Na-ue(> S) and these lead to a measurement of the total CIB 
intensity. 

Given an NEW model for each cluster in the survey, the 
magnification at each background pixel in the image plane 
map can be computed. The magnifications derived from the 
lensing model are applied to the measured source catalog to 
determine the source plane flux of each of the sources. The 
model is different for each cluster and the demagnification 
factor for each source depends on its position. Similarly, 
to compute the completeness of the survey, the source plane 
noise maps are divided by the same set of model magnifica- 
tions to create a background error map in which the value at 
each pixel is the effective measurement error at the source 
plane. These maps are multiplied by the source detection 
threshold value from Zemcov et al. (2007), and the resulting 
pixel values are summed to give the completeness as a func- 
tion of source flux, C(> S): the 50% completeness of this sur- 
vey, averaged over all fields, is ~ 5mJy in the source plane. 
The distribution of source plane fluxes is then corrected using 
C(> S) and summed into bins to yield A^meas(> S). 

The uncertainties in the derived A^meas(> S) depend on the 
measurement errors of the catalog, the uncertainty associated 
with the lensing models, and cosmic variance in the source 
sample. The most concise method of accounting for all of 
these sources of uncertainty is to use a Monte Carlo approach. 
The input parameters are varied over their nominal uncer- 
tainty range; these lead to a set of A^meas(> S) which span the 
allowed range of output values given the input uncertainties. 
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Determining the uncertainty due to cosmic variance is 
straight forward: in each simulation, a random realization of 
sources consistent with the counts is assigned uncorrected 
random positions before the lensing model is applied. Over 
many such simulations, variations in A'meas(> S) due to cos- 
mic variance should be measured to high accuracy. In each of 
the uncertainty simulations discussed below, cosmic variance 
is included. 

For each source, the IZemcov et al.l (l2007h catalog assigns a 
measurement error whose effect on A^meas(> S) is determined 
as follows. Each source's nominal flux is varied by an amount 
drawn from a Gaussian probability density function (PDF) 
with standard deviation equal to the measurement error of the 
source. This modified flux is then demagnified according to 
the model to determine its equivalent source plane value. 

Estimating the full uncertainty requires that the lens model 
parameters be varied on a per cluster basis. One of the least 
constrained parameters is the source redshift. The difficulty 
associated with determining unique counterparts of sub-mm 
sources at other wavelengths, whether in galaxy clusters or 
not, has rendered the precise redshifts for most of the catalog 
sources uncertain or unmeasured. As lensing magnifica tion is 
not a strong function of source redshift beyond z ~ 1 (iBlainI 
Il997h . the approach utilized in this analysis is to draw from 
a source redshift PDF based on the n{z) results presented in 
IXretxaga et al. (2007). In each calculation, the source plane 
for each cluster is assigned a redshift drawn at random from 
this distribution. Over many calculations, the full deflector- 
source redshift Zds parameter space is probed. The set of Zds 
could be varied on a per source basis; here, for computational 
efficiency, we choose to vary it per field. Over many realiza- 
tions this is an equivalent procedure. As these are well-studied 
clusters with precisely known redshifts, the cluster redshifts 
are never varied. A second parameter which is varied is the 
cluster's velocity dispersion, cr,,. These are typically available 
in the literature and are used here as the mean and standard de- 
viations of Gaussian distributions from which cr,, is randomly 
drawn; the scatter in the results of such realizations is a mea- 
surement on the effect of uncertainty in the lens model on 

A^meas(> S). 

Though each of these parameters could be varied alone or 
in any combination, all are varied simultaneously as we deter- 
mine the uncertainty in the measurement of A^meas(> S). This 
accurately reflects the lack of covariance between the different 
model parameters. In the results presented here, 10^ simula- 
tions were performed and propagated through to A^meas(> 5'). 
In each flux bin, the uncertainties in A^meas(> S) are deter- 
mined by calculating the set of A^(> S) which encompass 68% 
of the simulation results about the most likely value and tak- 
ing the extrema of the set. 

Many previous source count results have been derived from 
detailed lensing models of clusters. However, these models 
are observationally expensive and time consuming to con- 
struct and, unlike simple observables like velocity disper- 
sions, are typically not available in the public realm. Two 
clusters in our sample which do have publicly ava ilable lens- 
ing m odels written using the LENSTOOL c ode (Jul io et aP 

f 37l) ai-e Abell 1689 (Limousin et al.l l2007l) and Abell 2218 
fasdottir et al.ll2007b . To compare the statistical properties 
of the NEW model to these more detailed models, we perform 
simulations where backgrounds, populated using random re- 
alizations of the source counts in the way described above, are 
lensed with both the LENSTOOL and NEW models. The set of 
magnified source fluxes (of which only a very few are magni- 
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Model Magnification Linear Fit Coefficient 

Fig. i. — Histogram of linear fit coefficients arising from a comparison 
of tlie magnified fluxes derived from the NFW model used in this work and 
the LENSTOOL models which are publicly available for Abell 1689 and Abell 
2218. These distributions, derived from 500 simulations, are scattered around 
1.0 with a standai-d deviation of ±0.07 for Abell 1689 and ±0.04 for Abell 
2218. This shows that the NFW model can be used as an accurate and unbi- 
ased statistical description of a set of magnifications on many backgrounds in 
a given cluster 

fled by a factor greater than a few) can be fit to one another to 
yield a linear scaling between the resulting flux populations. 
Over many realizations, if the NFW model is capturing the 
statistical behavior of the lensing magnification as modeled 
by the LENSTOOL model, the scaling between them should be 
equal to I . As the NFW model can lead to very large mag- 
nifications, we set a threshold maximum magnification of a 
factor of 25; sources whose magnification factor is larger than 
this are set to this limit. Figure [T] shows the results of 500 
such simulations for both Abell 1689 and Abell 2218 models; 
the mean relation is indeed consistent with unity and the de- 
rived fit slopes have standard deviations of about 5% in both 
clusters. Uncertainties in the input velocity dispersions move 
the mean of the distributions within this standard deviation; 
presumably, a biased velocity dispersion measurement would 
also bias the lensing model though the bias would have to be 
> 100 km s"' to have much effect. This is evidence that, when 
used in Monte Carlo simulations, the NFW model accurately 
captures the bulk magnification properties of the more expen- 
sive models. 

2.3. Flux Bias Correction 

Though computation of A'meas(> S) and its errors does not 
require any a priori assumptions about the shape of Mriie(> S), 
correcting A^meas(> S) for the effects of flux boosting to deter- 
mine A^true(> S) does; the effects of flux boosting are estimated 
as follows. As a model parameterizing A^(> S) we use the bro- 
ken power law 



dN 



S 



— =No(- 
dS \So 



for 5 > So, 
for 5 < So. 



(1) 



A^(> S) is the cumulative integral of Equation ([T) from infinite 
flux to flux S. This broken power-law functional form for 
A^(> S) has a long pedigree in SCUBA blank field studies and, 
although not physically motivated, fits the current data above 
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the confusion limit very well (see, e.g. JCoppin et alj (l2006l) ). 
Since the total source plane area of this survey is not large, 
there are few intrinsically bright sources in it. To construct 
a model for source de-boosting, we therefore fix the model 
parameters which are best measured using sources w ith S > 
lOmJy , namely (3 and ^o, to the values found by Coppin et alJ 
( I2OO6I) . 5.1 and 9mJy, respectively, and only ever vary a and 

At a fixed a and A^o a set of sources whose flux distribu- 
tion is consistent with Equation ([T]) are assigned Gaussian dis- 
tributed random positions within a square grid. These maps 
are then lensed using the set of NEW cluster mass profiles. 
The calculation of flux boosting in SCUBA data is compli- 
cated by the fact that SCU BA is a double-differencing instru- 
ment ( Coppin et alj2 006). To account for this, each simulated 
source's flux is modified according to: 



(2) 



Here S is the flux in the simulated image plane map at po- 
sition (a, (5) where 'c' denotes the on source chop position 
and T and 'r' denote the two chopped positions, cr is a re- 
alization of noise drawn from a Gaussian PDE with standard 
deviation equal to that measured in the SCUBA map, and S' 
is the biased flux of the source. To reduce the complexity of 
these simulations, the SCUBA chop thr ow is set to 60'^ which 
is the median of the values used in the IZemcov et al.l (l2007h 
sample. The flux of each source in the simulated catalog is 
varied according to this algorithm, and the simulated source 
counts A^sim(> S) are formed by summing into bins in the usual 
way. As the simulated sources are boosted in the same way as 
the real sources and the simulated and measured catalogs are 
identically binned, A^sim(> S) can be directly compared with 
the boosted measured source counts A^mea,s(> S). 

Simulated, flux-boosted counts are calculated for a range 
of the parameters a and Nq; the outputs of such simulations, 
N^^''(> S), are then directly compared to A^meas(> S). Eor 
each pair of {a, No), 100 simulations are performed to mea- 
sure the mean and the scatter in A^^,^°(> S) arising from the 
variation in different noise realizations of Equation (|2]); the 
variance in the A^s"^°(> S) is typically 2 orders of magnitude 
smaller than the uncertainty in A'nieas(> S) and does not sig- 
nificantly contribute to the overall uncertainty. For each mean 
C*(>5)value, 



X = > , A^meas 



OS) 



-A?"'^°(> S) 

Sim ^ ' 



.2 

meas ' 



(3) 



where the sum is over flux bins. This is converted to a like- 
lihood function which allows us to assess the agreement be- 
tween the different model realizations and A^meas(> S). 

To correct A^meas(> S) for flux boosting, the model with 
(a, No) yielding the most probable fit to the data is determined, 
and for each flux bin the ratio of the boosted to input A^(> S) 
is calculated and used to correct the Nmeas(> S) to determine 

Mi-ueO S). 

3. RESULTS AND DISCUSSION 

The isotropic cumulative source counts A^true(> S) for fluxes 
below lOmJy inferred by our catalog are shown as filled dots 
connected by a solid line in Eigure [2] along with prior mea- 
surements of the density of 850 /im selected sources. Table |2] 
lists the measurements plotted in Eigure |2] for reference. Our 
results are consistent with prior data taken above SCUBA's 
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Fig. 2. — Cumulative source counts A'(> 5) derived from the lZemcov et alj 
f2007h catalog using the analysis described in the text (red points). The error 
bars account for the uncertainties associated with source flux measurement er- 
rors, gravitational lensing model uncertainties, flux boosting and cosmic vari- 
ance. The other points are from the fo l lowing sourc es: Small etal. l200j) 
Cowie et al. 1 2002), Webb et al. l2003h . lKnudsen et alJ <2006i) . ICopDin etalj 
1,200^, and . Knudsen et al. qOOS) . 



TABLE 2 

A'(> S) Measurements from this Work. 



S (mJy) 


N{> S) (deg-2) 


<JN+ (deg -) 


(7N- (deg ^) 


8 


259 


41 


45 


6 


541 


68 


74 


4 


974 


119 


109 


2 


1950 


340 


310 


1 


3730 


950 


720 


0.5 


7620 


2410 


1810 


0.25 


17300 


6700 


5000 



confusion limit, however they are lower than previously re- 
ported work for S < 1 mJy. 

Our catalog surveys a larger set of clusters than any previ- 
ous work, with a factor of ^ 2 more area at the source plane, 
and thus produces statistically tighter error bars. Importantly, 
the sources in this catalog are sampled from clusters which 
were not necessarily selected for their strong gravitational 
lensing. We believe our data are less biased than some pre- 
vious lensing surveys. Several of the targets - Abell 2218 
and MS045 1 .6-0305 are examples - are known to contain 
strongly lensed sources. Previous work is concentrated on 
these and similar fields to yield the largest number of sources 
per cluster with a median number of sources per cluster field 
of 4. In contrast, the catalog used here has a median number 
of two sources per cluster field. Another difference from pre- 
vious work is the conservative data ana lysis and point-source 
extraction used in IZemcov et al.l (l2007l) which led to the cat- 
alog used here. In that work, it was found that many sources 
from previous analyses drop significantly below the selection 
threshold. 

To check whether systematic artifacts dominate the differ- 
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1.0 10.0 
850 /xm Flux Density S (mJy) 

Fig. 3. — Compai'ison of counts measured near 850 and model pre- 
dictions. Tlie diamonds sliow tlie points meas ured in this work, tlie tri- 
angles are the counts measur ed by SHADES iCoppin et^ljl2006) and the 
dots are the counts measured by LESS i WeiB et al. 2009) where the source 
flux axis has been scaled by a factor (870/im/850/im)'' (amounting to 
a few percent) to account for the difference in the centroid wavelengths of 
SCUBA and LABOCA. The solid line shows the best fit of Equation Q] as 
described in the text. The dashed line shows the prediction using the model 
of Lagache et al. (20 031) , and th e dot-dashed line shows the prediction for the 
model of Valiante et aD f2009"). Both models require a break at ~ 1 mJy to 
obtain a total CIB consistent with the EZRAS result. The break in counts 
at Si 8mJy hinted in the LESS and SHADES work and presented here shows 
that the unverified second break in the counts at 1 mJy is not required to match 
the total CIB intensity. 

ence between this work and prior number counts, we have di- 
vided our catalog in several ways: into subsamples based on 
catalo gs presented in pr evious work (that of Small et al.ll2002l 
and iCowie et al.l 120021) . into targets with integration times 
longer and shorter than 15ks; and the analysis is repeated on 
these subsets. We find that the subsamples are consistent with 
the results presented for our full catalog and are lower than 
those from previous work. This is evidence that the subset of 
very deep or previously studied clusters have little statistical 
difference from the overall set. 

The 850 /im counts found here are plotted in Figure 
[3] along with the most preci se m easurements at brighter 
fluxes from SHADES (Coppin et al. 2006) using SCU BA and 
LESS using the LABOCA instrur nent at 870 ^ni ( WeiB et all 
20091). Model pred ictions from iLagache et all hoOSh and 
Valiante et al.l (l2009l) are also shown'. The shapes of the 



SCUBA and LABOCA points are very similar, though the 
LESS source density is a bit lower. The bright end of the 
lensed counts presented here are higher than either of the di- 
rect surveys, and we believe that our brightest few points may 
be biased high due to misidentification of a handful of clus- 
ter member sources as weakly lensed objects. Our survey 
is an unbiased count of lensed objects since clusters are not 
expected to be particularly aligned with background sources, 
but clearly the number density of cluster members is a bi- 



These 



models 



available 



ased estimator of the general isotropic source density. Incor- 
rectly counting cluster members as unit-magnification back- 
ground galaxies would inflate the brightest counts above the 
pure background-only count locus. 

The data support a single break in a power-law distribution, 
as in Equation ([TJ, with the break occurring at S'o = 7 or 8 
mJy. At the same time, it is clear that the full data sets are not 
consistent with each other given the quoted uncertainties. In 
fitting Equation (HJ to the data, we hold /3 fixed at 5.1, taken 
from Coppin et al., choose and then fit A^o and a to our 
five lowest flux data points. In this procedure, the choice of 
So essentially sets the level of the brightest counts with 8 mJy 
running through the SHADES points and « 7 mJy matching 
LESS. Reasonable choices of So have approximately a 1% 
effect on the inferred total flux of 850/im sources. Our best 
fit for = 8mJy is shown as a solid line in Figure |3] 

The total brightness of the sky contributed by all sources 
brighter than is 



NoSl 



l-(Si/So)^ 



13-2 



(4) 



online 



http : / /www, ias ■ u-psud. f r/ ir galaxies /model -php] and 
http : //www .physics . ubc . ca/^valiante/model| respectively. 



(where a = 2 would be treated as a special case). Values 
of vliiy) corresponding to our best-fit value of a = 1.86 and 
to the Icr error limits on a are listed in Table [3] in units 
of nW m"^ sr"' for several choices of the lower limit in 
flux. Si. A summary is plotted in Figure |4] The FIRAS ex- 
periment on COBE (Mat her et al.l 11994) measured the total 
extragalactic background in t he SCUB A band to be i n the 
range 3-5 x 10"'°W m"^ sr"' (iFixsen et a l. 1998; Pug et et al.1 
Il996l) . Given that only approximately half the total intensity 
of the CIB at 850/im comes from sources with S > 100 /iJy 
accessible to SCUBA lensing surveys, our constraints on the 
total brightness of the CIB from sources do not add to the pre- 
cision available from FIRAS for the total intensity of the CIB. 
At the same time, it is clear that the data do not require any 
break in the observed low-flux power law in order to be con- 
sistent with the FIRAS limits. Previous source number count 
determinations using SCUBA have the property that the ex- 
trapolation of A^(> S) over all fluxes overproduces this total 
CIB by a significant amount; in order to avoid this problem, 
either a second sharp break or a shallower faint count index 
would have been needed. 

Predictions given by the models of the infrared background 
source population shown in Figure [3] qualitatively reproduce 
the observed source density, steep slope, and break in the 
power law. However, they exhibit a break at ~ 1 mJy which 
follows from matching the evolution of the sources to the 
counts throughout the infrared regime. Altering simply the 
luminosity evolution terms in these models could match the 
faint counts but would not move the break in the power law to 
match the data. More thorough exploration of the combined 
effects of the luminosity and density evolution to understand 
disagreement between models and observations should be ad- 
dressed in future work. 

In addition to the counts determined here at 850 /im, source 
counts from the BLAST instrument working at shorter sub- 
mm wavelengths also show a flattening of the faint end of 
the curve (Patanchon et al. 2009). Furthermore, deep obser- 
vations of a lensing cluster with LABOCA independently 
suggest that the S^jOfj,m < 2mJy pait of th e A^(> S) curve 
has a flatter slope than previously measured (iJohansson et al.l 
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TABLE 3 
CIB Fluxes Varying a. 



a 


No 


ul all S 


5 > lO/iJy 


5 > 100/iJy 


S > 250fiJy 
(nWm-- sr') 




(deg^- mJy"') 


sr') 


sr') 


(nW sr' ) 


1.78 


98.0 


0.340 


0.269 


0.221 


0.194 


1.86 


92.0 


0.494 


0.309 


0.238 


0.203 


1.92 


86.8 


0.848 


0.364 


0.266 


0.222 



This Work. All S 
This Work. S > 0.1 mjy 
Coppin et 01. (2006), S > 2 mJy 
Smoil et 01. (2002). S > 0.15 mJy 
■ Fixsen et ol. (1998) 
Puget et ol. (1996) 




2.0 



1.5 



1.0 



0.5 



0.0 

12 3 4 5 6 7 

(x 10"'° W m"' sr"') 

Fig. 4. — Likelihoods of selected measurements of the extragalactic back- 
ground at 850 /im; those measured by integrating various estimates of the FI- 
RAS extragalactic background spectrum over the SCUBA band are shown as 
solid lines, while lower limits on the CIB from integrating the SCUBA source 
counts are shown as dashed. The red dashed line shows the limit on the CIB 
resulting when the counts presented here are integrated to S > 0. lOmJy; the 
PDF has a width which is the quadrature sum of the uncertainties on the 
best fitting a and Nq described in the text. The dashed black fine shows 
the limits which result when the best-fitting N( > S) model is extrapolated to 
5 = OmJy. The measurement of 'Copoin et al.' f2006), which integrates the 
number counts to S > 2 mJy, and a representative result from earlier SCUBA 
surveys of lensing clusters integrating the counts down to S > O.lSmJy are 
shown; because their power laws are steeper, these yield much larger intensi- 
ties if their best-fit models are extrapolated to very faint flux. 



I2OIOI) . These results corroborate our suggestion that the slope 
of the faint end of the sub-mm source counts has been overes- 
timated. 

The key to future measurements attempting to resolve the 
850 /im background into individual sources will be to survey 
large numbers of galaxy clusters to find rare, highly magnified 
image s of very dim sources. Experiments such as SCUBA- 
2 (.Holland et al., 2006) will allow mapping of > 100 clus- 
ter fields to the confusion limit of the JCMT in a matter of 
weeks of observing time. This can be compared to the perfor- 
mance expected from ALMA, which will have much better 
angular resolution than any single-dish instrument but a very 
small field of view. At 350 GHz, ALMA can generate a mo- 
saicked map of the central r = 2' of a cluster in about 100 



pointings, each requiring ~ 3 minutes to complete. In con- 
trast, SCUBA-2 can make confusion noise dominated maps 
of clusters far beyond this central region with an instrument 
noise of 0.5 mJy in 1 hr. Per cluster, this is a raw mapping 
speed advantage of a factor of 5 over ALMA on the central 
area; a survey strategy which would be maximally efficient 
would be to find bright sources with SCUBA-2 and follow up 
the central regions near the caustic lines with ALMA when 
interesting sources have been identified. 

Should a survey from a bolometric camera find bright 
sources in many clusters, it will be difficult to apply very de- 
tailed lens models to all of them, particularly ones found in 
upcoming millimeter SZ effect surveys which may not have 
ancillary data on which to rely. Statistical approaches like or 
more mature than the one adopted in this work will be useful 
for analyzing those data. More effort to create high fidelity 
lensing models - even if only accurate on average - based on 
simple observable quantities is certainly justified. 

That said, precise mass models of individual clusters which 
reduce the uncertainties on the source plane properties of 
highly magnified sources remain a very important component 
of these studies. Simple models are not very accurate when 
studying individual sources, so it is certainly worth expending 
the time and effort required to generate detailed lensing mod- 
els. Combined with such models, the data produced by very 
high angular resolution instruments like ALMA will surely 
allow us the best view of the faintest galaxies making up the 
sub-mm background. 
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